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ABSTRACT 



Since a nonlinear feedback control system may possess more than one 
type of forced oscillations, it is highly desirable to investigate the 
type of forced oscillations which can occur when the nonlinear restoring 
force function is of a specific form. 

This paper is to present a harmonic linearization method for finding 
the existence of forced oscillations and response curve characteristics 
of a nonlinear feedback control system by means of finding the restoring 
force function of the nonlinearity and using the harmonic balance and an 
iteration method for Investigating the conditions for one type of forced 
oscillations exhibited. 

The existence conditions for fundamental frequency, sub-harmonic of 
2nd order and 3rd order forced oscillations of a second order feedback 
control system are investigated; also the fundamental frequency forced 
oscillation for a higher order system and the jump resonance frequencies 
of response curve are investigated, and a general expression for the equi- 
valent gain of the nonlinearity has been developed. 

The author wishes to express his appreciation for the assistance and 
encouragement given by Dr. George Julius Thaler of the U. S. Naval Post- 
graduate School in this investigation. 
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CHAPTER I 



INTRODUCTION 

One of the major branches of servo mechanism analysis for nonlinear 
control systems is the study and methods of prevention of continuous oscil- 
lation. For a linear control system, the output may oscillate with con- 
tinuous or increasing amplitude with no input signal applied, in which case, 
the system is said to be unstable, or they may have oscillations which die 
away. In both cases this transient response is the characteristic of the 
system Itself and the conditions for divergent, continuous or decaying 
oscillation do not depend on the form or magnitude of the input signal. 

When a nonlinear element is present, this independence of the input 
signal r.o longer holds. The characteristic performance of the nonlinear 
control system will depend on both, i.e., forcing function and the 
characteristic of the nonlinear element. 

A nonlinear control system, under suitable conditions may exhibit 
steady oscillations in which the main component has a frequency which is 
dependent on the frequency of the forcing function. This type of oscilla- 
tion is called "Forced Oscillation". A forced oscillation for which the 
frequency is a fraction of the forcing function frequency is called "Sub- 
harmonic Forced Oscillation". A forced oscillation for which the frequency 
is multiple of the forcing frequency is called "Super Harmonic Forced Oscil 
lation". 

Since a nonlinear control system may possess more than one type of 
forced oscillations, it is highly desirable to investigate the type of 
forced oscillations which can occur when the nonlinear restoring force 
function is of a specific form. 

Some methods of investigation for nonlinear forced oscillation have 
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been discussed by J. C. West^ In his book "Analytical Techniques for 

Nonlinear Control System" (1960). A graphic method for investigating 

forced oscillation of "ON OFF" control system has been developed by Hamel 

2 

(French) and Tsypkin (USSR) . Some work for subharmonic oscillations of 

a typical restoring force function has been done by C. A. Ludeke and 

3 o ^ 

William Pong in their paper 1959^ and also by Ogata in his PhD thesis , 

June 1956. 

There are some other authors for investigation of the vibration of a 

5 6 7 

mechanical system, such as N. Minorsky , J, J. Stoker , and Y. H. Ku in 
their books of **Introduction to Nonlinear Mechanics” ;**Nonlinear Vibrations” 
and ”Analysis and Control of Nonlinear System” respectively. 

The purpose of this paper is to present a harmonic linearization 

g 

method for finding the existence of forced oscillations and response 

curve characteristics of a nonlinear feedback control system, by means of 

9 10 

finding the restoring force function of the nonlinearity ® , and using an 

iteration method for investigating the condition for one type of forced 
oscillations exhibited. A response curve and ”Jump Resonance” can be 
also investigated. 

The emphasis of this paper is placed on showing the general approach 
of the method, therefore, second order nonlinear feedback control systems 
are used as examples. However, a higher order system will be discussed in 
this paper also. 

By the harmonic linearization and iteration method, a general equa- 
tion for the equivalent gain of an odd function nonlinear element will be 
developed^ It is in terms of the amplitude and the frequency of input. 

There are six chapters in the main body of this paper. The first 
chapter and the last chapter are a general introduction and conclusion 
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respectively. Chapter two is a general description of oscillation of 
feedback control systems and an investigation of the restoring force func- 
tion of the nonlinear element. Chapter three investigates the existence of 
fundamental forced oscillations and the jump phenomena in the nonlinear 
feedback control system. Chapter four is a sub harmonic oscillation in- 
vestigation. Chapter five investigates the existence of forced oscillation 
of a higher order feedback control system. 
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CHAPTER II 



GENERAL DISCUSSION OF OSCILLATION IN FEEDBACK CONTROL SYSTEMS 
2-1 General Description: 

For either a linear or nonlinear feedback control system, most re* 
quirements for control system design specify a minimum time response, high 
accuracy, and high stability. It is clearly in order to have a minimum 
time response and high accuracy for a feedback control system, but first 
of all, the system should be stable. 

In general, two types of oscillations may be exhibited in a feedback 
control system, one is free oscillation, i.e., an autonomous system* This 
is a system which is left by itself with its outside source of energy sudden 
ly cut off or abruptly changed in amplitude. In such a system, free oscilla 
tions occur as the system tends to equilibrium after reaching a state of un- 
equilibrium. The oscillation frequency of this type depends only on the 
characteristic of system itself. Ordinarily, a free oscillation is a damp- 
ed oscillation; as the time increases the amplitude decreases and damps out 
in steady state. For an ideal case, if it is a system without damping, i.e. 
oC 9 0^ the oscillation will be continuous. 

The other type of oscillation is a forced oscillation, i.e. a nonauto- 
nomous system. This is a system that is acted upon by an outside source of 
energy, a force function. The force function may be a constant; a periodic 
function or any other function of time. 

The analysis and methods for preventing oscillations of a feedback 
control system are concerned directly or Indirectly with the characteristic 
of system, i.e., the differential equations of the system. If the system 
is linear a frequency response of the system can be directly solved from 

f 

the differential equation, furthermore, the principle of superposition can 
be applied too. On the other hand, if the system is nonlinear, either with 
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a nonlinear damping or nonlinear restoring force function, then the basic 
tools for the analysis of a linear system are no longer valid* 

2-2 Oscillation in Linear Systems: 

There are two types of methods for analysis and design of linear feed- 
back control systems; one is graphic method, the other is an analytical 
method. The best known of graphic methods consist of Bode diagram; Polar 
plot; Nichols diagram; and Root locus. If we know the transfer function of 
the system, any of the above methods can be used for solving the problem of 
stability. The analytical method is a mathematical analysis for solving the 
differential equation, and plot the curves of the solved equation. 

First consider a simple linear feedback control system, as shown in 
Fig. 2-1. The forward element has a transfer function G(s), and a feedback 












cfe;) 






^ /*f '> ) 

L 



Fig. 2-1 Block Diagram of a Linear System 
element transfer function H(s). The feedback signal is H(s)0^, (s) and it 
is assumed that this is subtracted from the input signal 0^ (s) to form 
the control signal e(s); thus: 

e(s) « iV(s) - H(s) S)^(s) (2-1) 

The relation between the control signal and reference signal in the 
Laplace transformation form is: 

- T V h(sjg(3) (2-2) 
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The solution of equation (2-2) for a particular Input can be obtained 
by Laplace and Heaviside methods. This involves the determination of the 
roots of equation: 

1 + H(s)G(s) = 0 (2-3) 



In order to split the transfer function into the sum of partial fractions 
of standard form for which time solutions are known. If equation (2-3) is 
an n th order polynorainal in s, there will be n roots say to such 
that (s- S^) is a factor of the polynominal. Special cases arise in which 
some of roots are identical, and if m roots of occurs, (s - S^)”' is a 
factor of the polynorainal. 

In general, the roots will be a complex number^ as: 



S s \ CU 

r r r 



(2-4) 



in which and are both real numbers, and since the coefficients of 

r r 

the polynominal of equation (2-3) are all necessarily real, then the com- 
plex roots should be in conjugate pairs, Thus: 



(S - + (S-^^ 



j ^ ) 

J 



will be a factor of the polynominal. 

In determining the transient response to a forcing function input, 
these various factors become the denominator of the partial fraction re- 
presentation of the transfer function equation (2-2). The term of K^/ (s *S^) 
gives a time solution containing expoential term: 

r + 3 CO ^)t 

The complete solution is the sum of all such terms. 

The complex conjugate pairs give rise to an oscillatory term: 

e ^ f ) 

If the system has convergent response, all of these terms in the time solu- 

tion must converge and hence all roots S, to S must have negative real 

1 n 

parts. If any of these roots has a zero real part, then a continuous 
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oscillation is produced. If any of the complex roots has a positive real 
part, then divergent oscillation will be the result. 

To determine the stability of a system by using the graphic method, it 
is not necessary to find the roots of the characteristic equation. If the 
equation is known algebrically, the Routh method can be applied and also 
by the Nyquist and Bode diagram method a graphical solution can be obtained, 
Nyquist criteron graphical analysis is the most useful method to solve 
the stability problem of control systems. Rewriting equation (2-3) in the 
form of: 

G(s)H(s) —1 (2-5) 

The roots of equation of (2-3) are particular values of s which satisfy 
equation (2-5). Thus if the equation of G(s)H(s) can be mapped for continu- 
ously variable values of s, then those values for which G(s)H(s) » -1 can 
be determined. 

Assume a particular roots of ^ + j can be represented by a 

particular value of G(s)H(s), i.e. 

G(^ ^ + J^cpH(.4^ + - R(^V, 0,,) ^ niU, , ^ ) (2-6) 

This equation represents a vector in I vs R plane, the magnitude of vector; 

r^2 2 

M « J R + I (2-7) 

and the direction by an angle; 

f « Tan"^l/R (2-8) 

For a constant the value of M will be a locus in I*R plane as 
CO varies. This process can be repeated for different values of<^ r to 
obtain a family of curves as shovm in Fig. 2-2. 

This steady state frequency response locus of the open loop system 
can be obtained experimentally from the system by measuring gain and phase 
for a sinusodial input of varying frequency. It can be shown that is the 
steady state frequency response locus passes to the right side of (-1, 0) 
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j 

I 



Fig, 2-2 Locus of an Open Loop Transfer 
Function G(s)H(s) 

point, this system will be stable, on the other hand, if it passes 
through or to the left side of (-1, 0) point, the system will be unstable. 
This criterion is due to Nyquist, A typical curve is shown in Fig, 2-3. 

In the extension to nonlinear control systems, or more complex multi- 
loop systems, the single loop criterion is sufficient for the majority of 
needs where G(s) can represent the over all transfer function of several 
loops provided that G(s) is itself stable. 




Fig. 2-3 Nyquist Criterion for Three 
Characteristic States 
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2-3 Forced Oscillations in a Linear System: 

Assume a second order linear system block diagram is shown in Fig. 

2-1 ; in which: 

and: 

H(s) * 1 (2-10) 

From equation (2-2), the differential equation of system becomes: 

e +<^0 K0 » K© (2-11) 

c c c r 

Let the forcing function: 

® ^ CostO t (2-12) 

Equation (2-11) becomes: 

0 + t/ 0 + K0 « C Cosr^t (2-13) 

c c c 

in which C is a constant of value of FK. 

The solution of equation (2-13) consists of the sum of the solution 
of the homogeneous equation, (i«e. the free oscillation of the system) and 
the solution of the non-homogeneous equatlon« Let the solution of homo- 
geneous equation: 

e^(t) * e”^ Cos6U^ ./ 1-/ t + C 2 stn OOnJ I'f" t) (2-14) 

In which; 

(On » JY 

^ s ’i- 

s Cos J' 

The complete solution of equation (2-13) will be: 

e (t) » e Coscv jY'f t + C, sinu Jf/t) ^ J (2-15) 

^ ^ nv ^ 

Equation (2-15) Is obtained by a superposition of the free oscillation 
and the forced oscillation which varies from the action of external force; 

The frequency of forced oscillation Is the same as that of the external 
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(2-16) 



force, the magnitude of the forced oscillation is given by; 



I ( / *V7 ' ^1. ; LK^ 

The value of phase shift 0 related to the external force Is given 

by: 

Cose - - (2-17) 

, )" t ry^ V )^ K * 

Sine - ^ — (2-18) 

^ -ij ' 

In the case of positive damping, i.e* U y 0^ it is clear from equation 
(2-15) that after a sufficiently long period of time, the free oscillation 
is damped out and only the forced oscillation would be observed* 

In the case of no damping, i.e. » 0, the phase shift 0 is seen 
from equation (2-17) and (2-18) to be zero for and 7L for In 

other words, the forced oscillation is in phase with the external force if 
the forced oscillation frequency is less than the free oscillation fre- 
quency, and is 180 degrees out of phase with the external force when u) is 
greater than dOy) . 

In the case of » 0, equation (2-15) becomes: 

0 (t) « (C_ Cos^' t 4* SintO t) 
c ' 1 n 2 ’^n 

+ C Cos(cOt 4- 0) (2-19) 

There will be a superposition oscillation of two frequencies, if the value 
of CO 4^ one of which is the natural frequency, and the other is the 
frequency of external force. In the case of u) , the free and the forced 
oscillations have the same frequency, and the solution of equation (2-13) 
will be found: 

C 

0 (t) * C- CoscX/t -f Sinat - ^,t Sint^jt (2-20) 

cl 2 

In this case, the component due to the external force is no longer periodic, 
it is oscillatory with an amplitude that increases linearly with time. At 
this condition, it is a resonance phenomenon. It is very important in de- 
sign to avoid this phenomenon in a system with a periodic force function* 
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In the case of a damped system, i<»e., > 0, from equation (2* 16) for 

the steady state, the amplitude of forced oscillation is always finiteo It 
is possible to give, by the use of dimensionless variable, a more general 
significance to the expression for the amplitude of forced oscillation, let 
the amplitude: 

C 

|H| - M( -- ) (2-21) 

Therefore equation (2-16) becomes: 

M - — r (2-22) 

h J f 

Where M is defined to be the magnification factor* The extreme values for 

2 2 2 

M are attained for iO * 0 and ( cv / cC\^ ) * 1 - 2 j-" * If 1 - 2 0, 

2 

there is a maximum for cU » 0; if 1 • 2y’ > 0, and there is a maxi- 
mum for =: /// ® minimum for co^ 0* For small values of damp- 

ing coefficient, the frequency which produced the maximum amplitude is very 
nearly the natural frequency of the system. Fig. 2-4 shows the forced oscil- 
lation response curves of a linear system as a function of with 

various values of . 

2-4 Oscillations in Nonlinear Feedback Control System: 

As mentioned in section 2-1 of this chapter, when a nonlinearity is 
present in a feedback control system, the properties of proportionality 
and superposition are not valid. In the linear system, these properties 
Involve the existence of a transfer function and of characteristic fre- 
quencies proper to the system (frequencies at which the system tends to 
oscillate with an amplitude depending on initial conditions); in the case 
of nonlinear systems, the amplitude, like frequency, can depend at one and 
the same time on both the initial conditions and the system itself. Some- 
times, as in the case of limit cycle both frequency and amplitude of input 
are characteristic of the system and independent of of Initial conditions. 
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In many cases the analysts and design of such systems using linear 
theory may produce an excellent result. This general procedure is a 
type of linearization procedure, i.e., the actual nonlinear system is re* 
placed by a linear system, which approximates it, and the analysis and de* 
sign techniques are applied to the linear equivalent. The most useful method 
for nonlinear linearization is that of the ’’Describing Function” method. It 
is actually a first harmonic approximation method. 

The method for linearization used by this paper is a multi-harmonic 
linearization and Iteration method; first assume a harmonic solution for 
the system differential equation and insert it in the differential equation, 
and compare the coefficients of the same order hantionic terms. A detailed 
procedure will be discussed in the next chapter. 

As previously explained, a nonlinear feedback control system may posses 
some type of oscillation. Which type will exist depends on the restoring 
force function of nonlinear element and the amplitude of forcing function. 
Hence, before investigating the existence of forced oscillations, the re- 
storing forcing function of the nonlinearity should be investigated. 

2*5 Characteristic of Nonlinear Elements: 

Nonlinear elements which may be seen in feedback control systems can be 
represented by a combination of four fundamental concepts of nonlinearity, 
which are: (1) Saturation; (2) Variable gain; (3) Dead zone; (4) Hysteresis. 
These basic notions of nonlinearity can, in combination with one another, 
result in almost all types of nonlinearity in practice, as shown In Fig. 2-5. 
For example an ideal relay nonlinearity can be seen as a saturation nonline- 
arity which linear part with a value of infinite slope. 

The operational characteristic of a nonlinear element which operates in 
a control system depends on the operating condition, just as a vacuum tube 
for which the operating characteristic depends on the operating point of its 
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(b) Variable Gain 
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characteristic curves. As a saturation nonlinearity element with a small 



forcing function input, it operates linearly. On the other hand, if the 
forcing function is very large and with a high frequency, the actual opera- 
ting characteristic is like an ideal relay. Hence, the actual operating 
characteristic of a nonlinear element may be said to be a function of the 
amplitude and frequency of the forcing function. 

2-6 Restoring Force Function Investigation of a Nonlinear Element: 

Consider the block diagram of Fig. 2-6, in which the output is a func- 
tion of input, usually this function in a control system is called the 



Nonlinear 



e 



e 



o 



Element 



Fig. 2-6 Basic Relation Between Input 
and Output of a Nonlinearty 



ti 



Restoring Force Function". 



Assume the input ist 



e^ ■ E(t) 



(2-23) 



The output will be a function of input, as: 



e^ « f(E) (2- 

9 

From the numerical analysis and curve fitting process, the function 
of f(E) can be expressed; 



(2-24) 




+ a E 
n 



+ b E 



n 



1/n 



(2-25) 
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Equation (2~25) can be in the general form: 



f(E) - a + a Ei + b E^^^ (2-25a) 

° l--'l ^ i'-“l 

for this paper, the nonlinear characteristic is considered S3rmmetrical to 
the original point, that means: 

f(-E) « -f(E) (2-26) 

It is an odd function, the coefficients of even order of equation (2-25) are 
zero, it becomes: 

3 5 

f(E) - a^ + a^E + a^E + a^E + 

+ b^E + + bjE^'^^ + (2-27) 

For the saturation and variable gain nonlinearity, J* C. West^ uses 
the restoring force function as a form: 

f(E) = a^E + a^E^ (2-28) 

in which is greater or equal to zero, and a^ is either greater, or less, 
or equal to zero corresponding to the case of a "Hard Spring", "Soft Spring", 
and "Linear Spring" saturation nonlinearity characteristic respectively. 

A set curves for all values of a^ is shown in Fig. 2-7. 




Fig. 2-7 Characteristic curves of Eq. (2-28) 
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In order Co generalize for the condition of saturation nonlinearity 
with dead zone, equation (2-28) can be written with a constant as: 

f(E) « -assign E + a^E + a^E^ (2-29) 

if the dead zone Is less than one, the characteristic curve can be shown 
in Fig. 2-8. 



riD 

ly/0 




Fig. 2-8 Characteristic Curves of Eq. (2-29). 

As discussed previously, when a saturation nonlinearity is operating 
with a larger amplitude than the saturation voltage and a high frequency, 
the operating characteristic is like an ideal relay. On the other hand, 
under some suitable conditions a relay nonlinearity control system may ex- 
hibit some sub-harmonic forced oscillations, as proved by the paper of A. M. 

12 

Hopkin and K. Ogata . This result is the same as a saturation nonlinearity. 

Actually, the operation of a relay is not a perfect discontinuous 
characteristic element, when the relay is operating from static to closed, 
there will be a little amount of time delay from starting to close con- 
tacts. Therefore, the operating characteristic of an ideal relay may be 
considered a saturation nonlinearity with a very large slope of the linear 
part, a sketch of it is shown in Fig. 2-9. 
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Fig. 2-9 Operating Characteristic of 
an Ideal Relay 

From the above discussion, and for a proper value of constants, the 
restoring forcing function can be expressed as: 

f(E) = a^E + a^E^ + a^E^ (2-30) 

Similarly, if we consider the dead zone the equation of (2-30) can 
be expressed: 

f(E) = -a^ SlgnE + a^E + + a^E^ (2-31) 

13 

M. J. Abzug investigates the restoring forcing function of an ideal 
relay nonlinearity by a fifth root or cube root of the input as; 

f(E) » bjE^'^^ (2-32) 

or; 

f(E) » b^E^'^^ (2-32) 

A more detailed discussion about frequency response by using the re- 
storing forcing function will be covered in the next two chapters. 
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CHAPTER III 



FUNDAMENTAL FREQUENCY FORCED OSCILLATION IN NONLINEAR CONTROL SYSTEMS 
3-1 General Description: 

Consider a system shown in Fig. 3-1, in which G^(s) is the transfer 
function of a controller element. G 2 (s) is the transfer function of a 
motor and gear element, G^ is a nonlinear element. H(s) denotes a transfer 
function of a feedback network. 

When a periodic forcing function is applied to the input, the output 
of the system may be or may not be a periodic function. If it is a peri- 
odic output and the frequency is the same as the frequency of the input 
forcing function, it is said that the system is operating in forced oscil- 
lation with the fundamental frequency. 



V ' 

y 







Fig. 3-1 Block Diagram of a Feedback 
Control System 

There are some methods to represent the forced oscillation in a non- 
linear system. If we keep an input forcing function with constant frequency 
and vary the amplitude of input, a response curve can be shown in the **In- 
put vs Error*’ plane, Fig. 3-2, in which there is a cut off amplitude; that 

means when the input is smaller than 0 , there is no forced oscillation. 

r-m 

The other representation is with the input amplitude constant vary 
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Fig, 3-2 Forced Oscillation Represented 

in the E vs 9 Plane 
r 

the frequency of input. It is shown in **Error vs Frequency** plane. Fig. 
3-2a is an example, it is actually a closed loop frequency response of the 
system, 

tr 

\ 

/ 



\ 

Fig. 3-2a Forced Oscillation Represented 
in the E vs cO Plane 

Either represented in E vs 9^ or E vs co plane, both of them can be 
represented in a phase plane (E vs E) . A typical phase plane of forced 
oscillation from the computer is shown in Fig. 3-3, in which it is represent- 
ed by an ellipse. The size of it is changed with the frequency. 

3-2 Basic Equation for Forced Oscillation of 2nd Order Nonlinear Feedback 
Control System: 

Consider a 2nd order feedback control system shown in Fig. 3-4: 





20 





I 








WV5 



C/„ o ) 









5C 



»* / 



Fig, 3-4 Block Diagram of a 2nd Order 
Feedback Control System 

The output equation in Laplace form: 






Kf -(E ) .-. 

s(S+c4 ) 



(3-1) 



in the differential form: 

e + 0 = Kf(E) 

c c 



(3-2) 



Where : 

0^(s) - 0^(s) - E(s) (3-3) 

Equation (3-2) becomes: 

E + ./: E + Kf(E) =0^ + ‘<^ 0^ (3-4) 

Assume the input forcing function: 

0^(t) »FCos(tt-t + 0) (3-5) 



in which 0 


is the phase difference between the 


error signal and forcing 


function of 


input. Inserting equation (3-5) in 


to equation (3-4): 






E + E + Kf(E) * —A Cos(ui t 


+ 0) — B Sin(iA^ t 


+ 0) (3-6) 


or: 


E + E + Kf (E) = C Cos 


(t^t+n + 9 + 0) 


(3-6a) 


Where : 


> 

H 

6 




(3-7) 




B a cl/ F 




(3-8) 




^ f.2 „2 

C * A 4“ B 


* a/u." C''" 


(3-9) 




0 ■ Tan ^B/A * 


Tan ^ c^/ to 


(3-10) 
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Equation (3-6) or (3-6a) is a basic differential equation of a 2nd order 
nonlinear feedback control system, In which f(E) Is the nonlinear restor- 
ing force function, and K is a constant of system. 

3-3 Forced Oscillation; Saturation Nonlinearity Without Damping: 

Consider a system without damping, l.e.^^ * 0, equation (3-6) becomes: 

E + Kf(E) « -A Cos(^<^t + 0) (3-11) 

Let the restoring force function: 

f(E) = a^E + a^E^ (2-28) 



Equation (3-11): 

E + K(aj^E + a^E^) « -A Cos(c:i/t + 0) (3-12) 

7 8 

By the harmonic iteration method * , the solution of equation (3-12) can 
be expressed as a series of odd order periodic harmonics: 

e(t) » Ej CoscJt + E^ Cos3tJt + E^ Cost<:0t ..... (3-13) 

in which E^, E^i are the amplitude of harmonics of e(t). 

If we consider only the fundamental frequency forced oscillation, we 
can assume one solution is: 



e(t) * Ej Cos ft) t (3-14) 

Where E^^ is the fundamental frequency amplitude of error signal, u> is the 
frequency of forcing function and forced oscillation; is to be determined. 
Since : 

E » - iO^E^ CosfDt (3-15) 

E » - u.' Ej SintOt (3-16) 

E^ » ^ E|^^(3 CosiOt + Cos3t'0t) (3-17) 

Inserting equation (3-14) into (3-17) In equation (3-12): 

2 3 3 1 3 

(- CO Ej + ^ 4 ^3^1 ) CostL^t + ^ a^E^ Cos3^^t (3-18) 

* -A (CosiUt Cos© + Sin^'Jt Sin6) 

neglecting the higher order harmonic terms and equating the coefficients 
of Cost^t and SinMyt, thus: 



23 



(3-19) 



(-(O^E, + Ka,E, + 7 Ka,El^) = - ACos6 
I 1 1 3 

ASln.vt Sine =0 (3-20) 

Equation (3-20), the term of ASln Ut Is not always equal to zero, the term 

of Sine should be zero. Hence the value of e will be either zero or «. 

Inserting equation (3-7) Into equation (3-19): 

Ka E T 3 Ka.E,3 

L! ) 1/2 (3-21) 



For e *! 0, and: 



a) = 



•Ka.El rj 
- ( 1 ^ 









(3-22) 



for 6 * Jt. 

As a check, in the case of F = 0, i.e, no forcing function input, equa- 
tions (3-21) and (3-22) give the free oscillation frequency of system: 

— 3 2 ^ 

cC - /K (a^ + ^ a^E (3-23) 

For the linear case, i.e., a^ « 0, the natural frequency : 

oOr, - (a^K)'^'" (3-24) 

It is exactly the same as we discussed in the linear system. 

Recall equations (3-21) and (3-22); the value of u) should be real, 
that means: 

^1 ®3^i^> ° 

for both the phase angle is either zero or n. And : 

Ej > F (3-26) 

for the phase angle is zero only. 

Equations (3-25) and (3-26) are the conditions for existence of the 
fundamental frequency forced oscillation. 

A free oscillation frequency response, i.e. F ® 0 can be sketched 
from equation (3-23), any value of force function other than zero also can 
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be sketched on both sides of the response curves of F * 0. It is to be 
noted, that the phase angle between the force function and the error signal 
is opposite, i«e,, 0 * n, when the response curves are to the left of curve 
for F = 0, and the phase angle between them is in phase, when the response 
curves are to the right of it. In other words, it is according to whether 
the frequency is less or greater than the frequency of free oscillation for 
that particular amplitude of error signal with a constant of force function 
input . 

A typical error vs frequency (E^>^ cu ) response curves and phase re- 
lation sketch for a different value of forcing function and different char- 
acteristic of nonlinearity are indicated schematically in Fig. 3-5. The 
response curve for free oscillation, corresponding F^O, is drawn as a 
dash line. 

In this respect the behavior of the nonlinear oscillation is the same as 
that of the linear oscillation. One sees that the response curve in the 
nonlinear cases could be thought of as arising from those for the linear case 
by bending the latter to the right for a hard spring saturation nonlinear- 
ity, and to the left for a soft spring saturation nonlinearity. 

A response curve in "Input vs Frequency" plane ( 6^ ) can be 

also obtained from relation of: 

e^(t) = e^(t) - e(t) (3-27) 

or: 

0 (t) « F Cos(a) t + 0) - E, Cosfoit (3-28) 

c 1 

In the case of 0 = 0, equation (3-28) becomes: 

0^(t) » (F - E^) Cosu;t (3-29) 

Hence : 

0^ = F - E^ (3-30) 
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or: 



E » F - 0 
I c 

Inserting equation (3-31) into equation (3-21): 

Ka (0^ -F) ♦ I Ka (0^ - F)^ 

fX' * ( s j 

For the case of G « it, equation (3-28) becomes: 

0 (t) « (F 4 EJ Cosa;t 
c 1 



Hence : 



1 



- (F + 



Inserting equation (3-34) into equation (3-23): 



K .^(e^ »3 (0, * F)^ 

0 



(3-31) 

(3-32) 

(3-33) 

(3-34) 

(3-35) 



From equations (3-32 and (3-35), if F « 0, the result is the same form 
as equation (3-23), that means the response curve for F » 0 is the same as 
plotted in the cO plane, only the difference is that the value of E^^ 
changes to 0^, The response curves for the values other than F ^ 0 are 
changed, and the phase relationship between the input force function and the 
output is reversed. 

There is a numerical example for plotting the response curve in 
and 0^ V plane respectively in section 3-7 of this chapter, and also 
experimental response curves are plotted. 

From the analysis of the above and the numerical example in section 3-7, 
it is shown that a saturation nonlinearity feedback control system can have 
fundamental frequency oscillations. This not only depends on the character- 
istic of the nonlinearity, but also depends on the amplitude of input force 
function. 

3-4 Forced Oscillations: Saturation Nonlinearity with Damping. 

As explained previously, if there is no damping present in the system, 
there is an oscillation error signal e(t) « EjCos(t(}t) either in phase or 
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out cf phase of 180 degrees with the forcing function of In the 

case where damping is present however, the forced oscillation displacement 
and the forcing input can be expected to be out of phase, just as in the 
corresponding of linear system. 

Recall the basic equation (3-6) and the restoring force function equa- 
tion (2-28): 

E +/E + Kf(E) * -A Cos(^'t + 0) - BSin( » t 4^ 0) (3-6) 

f(E) = a^E + a^E^ (2-28) 

Then equation (3-6) becomes: 

E 4 E 4 Ka^^E 4 * -ACos(ti/^t 4 0^ - BSin( t » 0) (3-36) 

Recall one solution of equation (3-36) is: 



e(t) s= Cos t 

Inserting equation (3-15) into equation (3-36): 

+ Ka^Ej + I Ka^Ej^^)Cos(^'t - ol u) 'E ^Sin u) t + ^ Ka2Ej^Cos3 



(3-15) 



t 



= (-ACOS0 - BSinG) Cosc^ t + (ASin6 -B Cosej .Sin t (3-37) 
Neglect the higher order terms of harmonic and equates the coefficient 



of Sin t and Cos<^ut respectively, hence: 

+ Ka,E, + 7 Ka-E,^) s ACos6 - BStn9 
1 11 4 3 1 

* A Sin0 - B Cos0 

Squaring equations (3-38) and (3-39), and adding: 

(- ti}''E, + Ka.E, + 7 Ka.E,^)^ ^ 

1 11431 1 



or: 



3 3222^ ^ 2 ' 2 

(-a.-Ej + Ka^E^ + ^ Ka^E^-") a) F t-' '■ F 



also: 



- F^) - A,"[2Ej^(Ka^ + | Ka^E^^) -.^■"(E^^ F^) 

3 

+ E^^ (Kaj^ + 4 Ka^Ej^)^ = 0 



(3-38) 

(3-39) 

(3-40) 

(3-41) 



(3-42) 
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As a check of equations (3-38), (3-39) and (3-41); first consider 
the response curve in the high frequency and the low frequency ranges, 
the phase relation between the error signal or control signal and the in- 
put forcing function is either zero or 180 degrees. In this case, the 
sine terms of equation (3-38) and (3-39) are zero, the response curve is 
the same as equation (3-20). That means for a damped system if the phase 
angle between the input and the output is either zero or 180 degrees, the 
frequency response is exactly the same as the system with no damping. 

On the other hand, if the value of damping is very small, the re- 
sponse curve will be very close to the sketch of those of Fig. 3-5. The 
only difference when very small damping is present is that the response 
curves are rounded off in the vicinity of curve for F » 0. A sketch of 
them is shown in Fig. 3-6. 

For the phase investigation; recall equation (3-39) and rewrite in 
this form: 



Sine - Cos0 sc>iE^/F 



(3-43) 



or : 






-1 



Sin(e - Tan ) = p ( (3-44) 

Hence, the phase angle between the error signal and the forcing function 



of input; 



0 » Sin 



•1 



+ Tan 



•1 



(3-45) 



F (f0' + ^) 1/2 a) 

Examine equation (3-45) in the case of » 0, 0 is either zero or Jt, 
when damping is present the value of 0 is a function of cO and if the 
values of and F are fixed. It should be noted, there are two values of 
phase shift for one value of error signal; one is for the low frequency 
( a. < ) and the other is for the high frequency one is from 

equation (3-45) and the other is from the 180 degrees minus the value of 
equation (3-45). 
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For the existence conditions for forced oscillation when damping is 



present, recall the equation (3-42), and put: 



3 2 

Ka, 4- 7 - Ka.E, » X 

1 4 3 1 



and: 



(Ej - -(2E^^X (Ej - 

Solve equation (3-47): 



2 2 

+ Ej X » 0 



(3-46) 



(3-47) 



U> ,j{ 1 'J U\br h (3-48) 

It is to be noted that the conditions for forced oscillations in an 
undamped system are still valid for a damped system within the low and 
high frequencies, in other words, it is valid for the condition of phase 
angle either zero or n. 

From equation (3-48) where forced oscillation exists, the value of 
should be real, hence; the condition for equation (3-48) to be real: 



and; 



or: 



2 2 
- F > 0 



2 .^2 2 

2E^X - - F ) ^ 0 



X . 






J. (E-^ - F^) > 0 



(3-49) 



(3-50 



(3-51) 



and; 



or: 



[2e,^X - - e') j ' - 4E^V(Ej^ - F^) > 0 



2„2,_ 2 _2. 



4 Ej^X^F^ + (Ej^ - F^) > ^Ej^X^ (Ej^ - F^) 



(3-52) 



(3-53) 



Equation of (3-49); (3-51) and (3-53) are the conditions for the 
forced oscillation of a damped system. 
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3-5 Jump Resonance: 

Recall a sketch of frequency response curve of a soft saturation non- 
linearity shown in Fig. 3-7. ^Then a experimental frequency response is 
being measured for a system, this measurement can be done by an analog 
computer to measure the amplitude of the input and the output, at the same 
time measure the phase relation between them. 




Fig. 3-7 The Region of Jump Resonance 
It is found that under certain conditions of amplitude and frequency 
an infinitesmal change of either frequency or amplitude of the input signal 
causes a large and discontinuous jump in amplitude of the output, at the 
same time, a discontinuous phase jump occurs. Fig. 3-8 shows a set of 
amplitude jump and phase jump of a typical system from an Electronic Analog 
Computer experimental work. 

Fig. 3-7 is shown as a typical system with a soft saturation non- 
linearity and with a fixed amplitude of the input, which is larger than the 
value of saturation of nonlinearity. If we are going to measure the frequency 
response, being with a low frequency, then increase the frequency, given a 
response from A to B. At point B with a frequency u)f , an infinitesimal 
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increase in frequency causes the ai^iplltude to jump upwards from point B to C, 
and with a phase shift from below 90 degrees to a value above 90 degrees# 
Further increases in the frequency cause the output to follow the gradual 
smooth curve to D. Now, if the input frequency is reduced from higher 
values, the output will follow the smooth curve DC. When the frequency is 
reduced to the point E, if the frequency is infinitesimal reduced, the out- 
put jumps down from point E to point F. 

It is found that either from the theory analysis or from the experi- 
mental work, a hard spring saturation nonlinearity can also produce the 
jump resonance, only the difference between them is the jump in the oppo- 
site direction from each other. A typical system for a hard spring satura- 
tion nonlinearity will be shown in the experimental section# 

From the above discussion, the portion of the cuirve from point B to 
point E cannot be obtained experimentally, usually this portion of the 
curve corresponds to an unstable operating condition# 

3-6 Lower and Upper Jump Frequency of the Jump Resonances 

As previously discussed, the jump frequency of 'O. is called ’’Upper 
Jump Frequency” and the jump frequency of oJ,^ is called ’’Lower Jump Fre- 
quency”. Recall the equation (3-41) and the sketch of the response curve 
shown in Fig. 3-7: 

(-(O' Ej + (3-41) 

The equation of the locus of the vertical tangents can be found by 
differentiating equation (3-41) implicitly with respect to and setting 
d (O equal to zero, the result is: 

(-a'-'- + Kflj + Ka2Ej^)(-^./+ Ka^ + Ka^E^^) (3-54) 

Check the equation of (3-54), if the damping is very small, the term 
Uj' U j ' can be neglected, that leads to a pair of equations: 

u.-- K(a^ + I* fl3Ej^) = 0 (3-55) 
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U K(a^ + -I- a3Ej^) “ 0 (3-56) 

The curves corresponding to these equations are shovm in Fig. 3-9. 

It is to be noted, the equation of (3-55) is the same as the response 
curve equation (3-23) for the free oscillation of an undamped system. The 
equation of (3-56) is the locus of the vertical tangents of curve of Fig. 
3-5. 

In the case when damping is present, the locus of vertical tangents 
equation (3-54) will be directed by the equations of (3-55) and (3-56), 
in particular there should be one branch near the response curve for free 
oscillation, and another near the points where the curves for the undamped 
system. In other words, the curve must appear as they are shown in Fig. 
3-10. 

Rearrange the equation (3-54), and solve for tC , which leads: 




2m, J T 5.' 






■ it\^ if 



■(. n 



(3-57) 



or: 












58) 



In which « = i,'Kaj> is the natural frequency of the linear system. If 
we know the constants of system, and the amplitude of fundamental fre- 
quency forced oscillation, the upper and the lower jump frequency can be 
calculated from either one of the equations (3-54), (3-57), and (3-58). 



3-7 Analog Computer Analysis for a Second Order System with Soft Saturation 
Nonlinearity : 

Consider a system with a block diagram shown in Fig. 3-11. 
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Fig. 3-11 Block Diagram of a 2nd Order System 
with a Soft Saturation Nonlinearity 

In which the restoring force function of nonlinearity is: 

f(E) = 2 e - 0.032E^ (3-59) 

The cliaracterist ic cuxrve of restoring force function is drawn in Fig. 3-12 

where : 

E * 1.00 1.50 2.00 2.50 3.00 3.50 4.00 5.00 6.00 

f(E)« 1.97 2.89 3.75 4.50 5.15 5.65 5.95 6.00 5.10 

Before the experimental analysis consider the theoretical calculations 
and curves. First consider the response curve on the Bi plane for the 
system without damping, and with an input function: 

0^(t) s 3 Cos(o;t + 6) (3-60) 

Recall the equations of (3-21); (3-22) and (3-23) and leads: 

= (4 - 0.048Ej^)^^^ (3-61) 

for free oscillation and: 




4E - 0.048e/ 1/2 

i. = (-* r 7-^3 ) (3-62) 

for the phase angle 9 * n, and: 



4E - 0.048e/ 1/2 

.'L = ( (3-63) 

E - 3 

for the phase angle 0*0. 
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The calculations and curves of equation (3-61) to (3-63) are shown 



in Table 3-1 and Fig. 3-13 respectively. 



Table 3-1; 


Calculation for 


response curves on the 


plane : 


Amplitude of 


Frequency for 


Frequency for 


Frequency for 


Error Signal 


Free Oscill. 


0 ss n 


0 a: 0 


( ) 


( - ) 


( ) 


( ) 


0.00 


2.00 


/ 


/ 


0.50 


1.996 


0.75 


/ 


1.00 


1.99 


0.99 


/ 


2.00 


1.95 


1.23 


/ 


3.00 


1.89 


1.37 


/ 


3.50 


l.f-4 


1.35 


4.88 


4.00 


1.80 


1.36 


3.60 


5.00 


1.68 


1.32 


2.65 


6.00 


1.53 


1.23 


2.16 


7.00 


1.29 


1.08 


1.73 


8.00 


1.00 


0.80 


1.22 


9.00 


0.35 


0.28 


0.42 



For the calculation of response curves and phase shift, when damping is 



present, recall equation (3-42), and leads: 

U.i-(r[0.25^- — 2 — 2" ■ 0-f48Bj )} + 

E - 3 



1 2 2 
A, (4 .0.048E:) s 0 

(3-64) 



for the phase shift, recall equation (3-45) 

-1 



0 IS Sin 






1/2 + Tan 



1 



(3-45) 



F( (0-^+ ) «. 

It is to be noted, from equation (3-45), and the response curve for an un- 
damped system shown on the e, a d plane, when the frequency is larger than 
the frequency of free oscillation with a particular amplitude of input, 
the phase angle 6 is zero and is rt, when the frequency is below the free 
oscillation frequency. Hence; the phase calculation from equation (3-45) 
should be modified as: 



0.25E, 



© - 180 



(Sin 



-1 



0.25 



3( >0+ 0.25^)^^^ 



+ Tan 



(3-65) 
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for c 



, and : 



e = Sin 



-1 



0.25E 



1 



+ Tan 



-1 



0.25 



2 t/'* 

3 ( a + 0.25 



(3-66) 



for Cl . 

If the response curve is presented on the cC| plane, the condition 
of equation (3-65) and (3-66) should be reversed. 

The calculation for the system with damping of 0,25 and the curves are 
shown in Table 3-2 and Fig, 3-14 respectively. 

The jump frequencies of system can be calculated from equation (3-57) 
and also can be seen from the response curve. It is shown schematically 
in the diagram of Fig. 3-14, 



Table 3-2: Calculation for Response Curves (with a damping of 0,25) 



Amplitude of 


Frequency 


Phase for 


Frequency 


Phase for 


Error Signal 


for u. 




for^^^^V't 




( c-/ ) 


( ^ ) 


( ) 


( cv ) 


( 0 ) 


1.00 


1.00 


164.00 


/ 


/ 


2.00 


1.24 


161.00 


/ 


/ 


3.00 


1.35 


158.50 


/ 


/ 


3.50 


1.36 


156.50 


4.83 


3.00 


4.00 


1.37 


155.00 


3.58 


7.00 


4.50 


1.40 


153.00 


2.80 


12.00 


5.00 


1.36 


151.50 


2.50 


15.50 


6.00 


1.24 


144.50 


2.05 


21.00 


7.00 


1.09 


135.00 


1.64 


28.50 


8.00 


0.87 


121.50 


1.16 


45.50 


8.70 


0.65 


90.00 


0.65 


90.00 



For the experimental analysis, the first investigating is the restor- 
ing forcing function. It is to be approximated as a two slope lines, one 
slope is 2, and the other is zero. The circuit setting in the analog com- 



puter is shown below, and also all components to be used are schematically 



in the circuit of Fig. 3-15, 

The characteristic curve for the restoring forcing function from the 
computer is shown in Fig. 3-16; (a) is shown the waveform of input (e^) 
and output (^^) , (h) is shown the characteristic of the output to the input. 
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Fig. 3-15 Operational Diagram of a Function 



Generator for f(E) * 2E 



0.032E' 



The operational diagram of analog computer setting up for system of 
block diagram in Fig. 3-11 is shown in Fig. 3-17, in which; 



e(s) = 






~ (s)j 



^(s) = 



) 



r / 



= 



^ ( 0 - 



(S) 

(s) 



From equations (3-67) to (3-69) system scaling is shown below: 



(3-67) 

(3-68) 

(3-69) 



where : 




= 1.00; 1.00; 


1,00; 


“ 2.00 






SEE 


®1 ^fl^^l “ 




« 1.00; 


Rj^ - 1.00; 


Rj = 1.00 




= 


®2 ^fl'^*^2 “ 


*2 


= 1.00; 


R2 = 1.00 




^3 




^f 2 ^fl/a^ = 4.00 




- 0 . 50 ; 


R^2 * 1.00 






■C 


*3^f2/V3 “ 


®3 


= 1.00; 


R3 * 0 . 50 ; 






- 


85/ R5CJ2 = 2.00 




= 1.00; 


Rj * 0.50 


Cj^2 = 1*00 


^^6 


SB 


a^R^ /R, = 1.00 
D t 0 


^6 


= 1.00; 


R, - 1.00 
6 


Rj2 ® 1.00 



in which R in megohm; C in ^ f 
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The results fiom the anelog computer experimental work Is condensed 
in Table 3-3; and the response curve v^ith phase shift represented on i, 
plane is shown in Fig. 3-18 respectively. Also the Jump frequencies are 
indicated schematically in these diagrams. 

A theoretically calculated response curve in /f, plane for comparing 
is shown with the experimental response curve of Fig. 3-18. Comparing the 
results from the computer experimental work with the results from theoretical 
calculation, it is shown that the results are very close within high fre- 
quency and low frequency range. There is a resonance phenomena between the 
jump frequencies. 
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Table 3-3 Data for Response Curves from Analog Computer for System 

* 0.25; K « 2; f(E) = 2E - 0.032 e3; £}.(t) - 3.0Cos(^t f 0) 



(a) Frequency Increasing; 



Frequency of 


Amplitude of 


Phase of 


Amplitude 


of 


Phase i 


Input 


Error Signal 


Error 


Output 




Outpui 


(O) ) 


( ^/) 


( ) 


( ) 




( <9 ) 


0.500 


0.38 


180.00 


3.20 




0.00 


0.628 


0.50 


180.00 


3.40 




0.00 


0.755 


0.75 


175.00 


3.60 




5.00 


0.880 


1.25 


168.00 


4.00 




9.00 


1.000 


1.50 


163.00 


4.50 




12.00 


1.130 


2.15 


158.50 


5.00 




24.00 


1.260 


3.00 


150.00 


5.40 




35.00 


1.260<+) 


9.00 


/ 


6.50 




/ 


1.380 


8.25 


25.00 


5.50 




140.00 


1.510 


7.30 


14.00 


4.50 




155.00 


1.640 


7.00 


5.00 


3.80 




159.00 


1.760 


6.50 


0.00 


3.40 




165.00 


1.890 


6.20 


0.00 


3.00 




172.00 


2.010 


6.00 


0.00 


2.60 




180.00 


2.140 


5.50 


0.00 


2.20 




180.00 


2.260 


5.20 


0.00 


2.00 




180.00 


2.390 


5.00 


0.00 


1.60 




180.00 


2.510 


4.80 


0.00 


1.40 




180.00 


2.830 


4.40 


0.00 


1.20 




180.00 


3. 140 


4.00 


0.00 


1.00 




180.00 


3.460 


3.80 


0.00 


0.80 




180.00 


3.770 


3.70 


0.00 


0.60 




180.00 


4.090 


3.60 


0.00 


0.50 




180.00 


(b) Frequency Decreasing: 










Frequency of 


Aniplittide of 


phase of 


Amplitude 


of 


Phase o 


Input 


Error Signal 


Error 


Output 




Output 


(cO ) 


( 


( - ) 


( Cl) 




( ) 


2.5’0 


4.75 


0.00 


1.50 




180.00 


2.200 


5.20 


0.00 


2.00 




180.00 


1.890 


6.20 


0.00 


3.00 




180.00 


1.760 


6.80 


0.00 


3.20 




173.00 


1.640 


7.10 


0.00 


3.80 




168.00 


1.510 


7.58 


0.00 


4.50 




162.00 


1.380 


8.30 


12,00 


5.50 




156.00 


1.260 


9.20 


23.00 


6.80 




140.00 


1.130 


10.00 


34.00 


8.00 




135.00 


1.000 


11.50 


65.00 


9.50 




120.00 


0.880 


11.00 


/ 


8.70 




/ 


0.880(-) 


0.75 


175.00 


4.00 




0.00 



48 






’> C 






’J 



-{ 

o* 





o 


















49 



The resonance phenomena in Fig. 3-' 18 is actually similar to the re- 
sonance in the linear system. It depends on the value of damping in 
the system. Another experimental result for the same system and same 
value of input, only changed the damping from 0.25 to 0.5, and is shown 
in Table 3-4 and Fig. 3-19 respectively. In this case the resonance pheno- 
mena is much lower than the case of a damping 0.25, 

It is to be noted, if the damping of a system is Increased the range 
of jump frequencies is decreased. For this reason, when the damping in- 
creased to limit value, at which no jump phenomena existed, it means that 
the system becomes a linear system. On the other hand, if the system is an 
undamped or a very low damping system, it will be very oscillatory or a 
limit cycle exists. 

The response curve represented on ^’Output vs Frequency** ) 

plane can be calculated from equations (3-32) and (3-35). An experimental 
response curve shown on the ^ oC plane from the system with a damping of 
0.25 is shown in Fig. 3-20. 

3-8 Analog Computer Analysis for a Second Order System with Hard Saturation 
Nonlinearity: 

Consider a system with a block diagram shown in Fig. 3-21: in which 



A 



1 








\ 






Fig. 3-21 Block diagram of a Second Order System 
with a Hard Saturation Nonlinearity 

assume the restoring forcing function of the nonlinearity is: 

f(R) = 0.425E + O.OOSE^ (3-70) 
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Table 3-^4. Data for Response Curve from Analog Computer for System 

= 0.50; K = 2; f(E) • 2E-0.032 e 3; 6,(t) * 3.0Cos(«.' t+o) 

(a) Frequency Increasing: 



Frequency of 


Amplitude of 


Fhase of 


Amplitude of 


Phase of 


Input 


Error Signal 


Error 


Output 


Output 


( ^ ) 


( i-, ) 


( ) 


( 4, ) 


( ^ ) 


0.628 


0.50 


173.00 


3.40 


7.30 


0.755 


0.75 


165.00 


3.50 


12.00 


0.880 


1.00 


155.00 


4.00 


16.50 


1.000 


1.50 


140.00 


4.75 


24.00 


1.130 


1.75 


125.00 


6.00 


30.00 


1.260 


8.30 


91.00 


8.00 


85.00 


1.380 


8.10 


47.00 


6.25 


120.00 


1.510 


7.50 


36.00 


5.25 


135.00 


1.640 


7.00 


25.00 


4.25 


148.00 


1.760 


6.50 


20.00 


3.50 


160.00 


1.890 


6.00 


12.00 


2.75 


172.00 


2.010 


5.50 


5.00 


/ 


/ 


2.140 


5.00 


0.00 


/ 


/ 


2.260 


4.50 


0.00 , 


1.75 


180.00 


2.390 


4.25 


0.00 


/ 


/ 


2.510 


4.00 


0.00 


1.25 


180.00 


2.810 


3.75 


0.00 


1.00 


180.00 


3.140 


3.50 


0.00 


0.75 


180.00 


3.460 


3.35 


0.00 


0.50 


180.00 


3.770 


3.25 


0.00 


0.30 


180.00 


4.090 


3.10 


0.00 


/ 


/ 


(b) Frequency Decreasing 








Frequency of 


Amplitude of 


phase of 


Amplitude of 


Phase of 


Inpul: 


Error Signal 


Error 


Output 


Output 


( a) ) 


( ) 


( ) 


( <5. ) 


( 6 ) 


2.510 


4.30 


0.00 


1.25 


180.00 


2.200 


5.00 


0.00 


1.75 


180.00 


1.890 


5.75 


0.00 


2.75 


172.00 


1.760 


6.25 


7.50 


3.50 


160.00 


1.640 


7.00 


15.00 


4. 10 


148.00 


1.510 


7.75 


27.00 


5.00 


132.00 


1.380 


8.25 


32.00 


6.50 


120.00 


1.260 


8.50 


60.00 


7.50 


97.00 


1.130 


7.00 


80.00 


8.20 


40.00 


1.000 


2.00 


130.00 


4.75 


19.50 


0.880 


1.25 


145.00 


4.00 


0.00 
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Fig. 3-21' Chc.ract3rlotic Cui'^/e of 
torin^' fviictlon : f v- ) — 0 .^!- 25 f ♦■ 0 . 0032 ^ 



The characteristic curve of equation (3-70) is shown in Fig. 3-2la, where: 

E « 1.00; 2*00; 3.00; 4.00; 5.00; 6.00; 7.00; 8.00. 

f(E) = 0.43; 0.89; 1.42; 2.02; 2.75; 3.63; 4.69; 5.96. 

A theoretical analysis can be calculated from equations (3-70) and 
(3-21) to (3-23), as we did in section 3-7. Only the analog computer anal- 
ysis is considered in this section. 

For the simulation of the restoring forcing function of nonlinearity 
from the analog computer, it is approximated by two straight lines, one 
with a slope of 0.425, and the other with a slope of 2.00. The circuit set- 
ting in the computer is shown in Fig. 3-22, and all components to be used 
are schematically in the diagram. 



^ - 



r, ^ ' 

.1 h I I- 






/ c 







Fig. 3-22 Operational diagram of the Function Generator 
for f(E) = 0.425E + O.OOSE^ 

The characteristic curve of restoring forcing function from the com- 
puter is shown in Fig. 3-23. 

The operational diagram of analog computer set up for the whole 
system is mostly the same with the diagram of Fig. 3-17, the only dif- 
ference is the function generator circuit. 
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The results from the analog computer are condensed in Table 3*5, 
and the response curves with phase shift shown on the ^ 

plane are shown in Fig, 3-24 and 3-25 respectively. In which the input 
forcing function is 0 (t) =3,50 Cos(^v t + 0). 
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Table 3-5. Data for Respoi se Curves from Computer for a System with a 
Fard Saturation Nonlinearity, In which: ^ 0.25; K » 2; 

f(E) = 0.425E 4 0.005E3; Q^(t) = 3.5 Cos(^ t 4 0). 

(a) Frequency Increasing: 



Frequency of 


Amplitude of 


Phase of 


Amplitude of 


Phase of 


Input 


Error Signal 


Error 


Output 


Output 


( (O) 


( ) 


( ) 


( tPc ) 


( (9 ) 


0.500 


1.75 


180.00 


5.50 


0.00 


0.628 


3.50 


180.00 


6.50 


4.00 


0.754 


6.50 


175.00 


9.00 


7.50 


0.880 


8.00 


168.00 


10.50 


10.00 


1.000 


9.00 


164.00 


12.00 


11.50 


1.130 


10.50 


158.50 


13.00 


13.00 


1.260 


12.25 


152. r>0 


14.50 


15.00 


1.380 


14.00 


145.50 


16.30 


17.00 


1.510 


16.25 


141.00 


18.75 


20.00 


1.640 


19.00 


134.00 


21.30 


24.00 


1.760 


22.25 


124.00 


22.00 


37.00 


1.890 


26.00 


102.00 


25.00 


45.00 


1.950 


26.50 


/ 


25.50 


55.00 


2.010 


26.00 


24.00 


25.00 


90.00 


2.140 


4.50 


20.00 


0.86 


150.00 


2.260 


4.35 


15.00 


0.65 


165.00 


2.510 


4.25 


7.50 


0.50 


180.00 


2.830 


4.15 


0.00 


0.40 


180.00 


3.140 


4.00 


0.00 


n. 30 


180 no 


^.770 


3.80 


0.00 


/ 


/ 


(b) Frequency Decreasing: 








Frequency of 


Amplitude of 


phase of 


Amplitude of 


Phase of 


Input 


Error Signal 


Error 


Output 


Output 


(.CO ) 


( 


( ) 


( ) 


( f: ) 


2.510 


4.25 


5.00 


0,50 


180.00 


2.400 


4.35 


■'0.5'^ 


0.63 


180.00 


2.270 


4.50 


13.50 


0.68 


172.50 


2.140 


4.60 


15.00 


0.75 


170.00 


1.890 


4.85 


21.00 


1.13 


157.50 


1.760 


5.00 


23.50 


1.25 


155.00 


1.640 


5.35 


25.00 


1.60 


152.00 


1.510 


5.75 


26.00 


2. 13 


150.00 


1.380 


6.75 


27.50 


16.30 


18.50 


1.320 


12.25 


/ 


/ 


/ 


1.260 


11.70 


150.00 


14.50 


15.50 


1.130 


10.25 


155.00 


13.00 


13.50 


1.000 


9.00 


160.00 


11.70 


12.00 


0.880 


7.50 


172,00 


10.50 


10.00 
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CH/'lrTER IV 



SUB“HAP>IONIC OSCILLATION IN NONLINEAR FEEDBACK CONTROL SYSTEM 
4-1 General Description. 

Up to the preceding chapter, there is only discussed the forced 
oscillations for which the frequency is the same as that of the external 
force. A nonlinear control system, under suitable conditions may have a 
steady state oscillation in which the main component has a frequency which 
is a fraction 1/n (n is a real and positive integer) of the forcing func- 
tion frequency. These oscillations are called *’Sub-Harmonic Oscillation'*. 

Recall a linear system (Section 2-3), there are two oscillation terms 
in the transient response, if the frequency of the free oscillation is 
/n (n is an integer), then a force function of frequency can excite 

the free oscillation in addition to the forced oscillation of frequency c . 

Wt/2 

But as t increases, e damps the free oscillation until only the 

steady state oscillation remains. Hence there is no sub-harmonic or super- 
harmonic oscillation existing in a linear control system at steady state. 

In the case of nonlinear control system, the steady state response 
may not have the same frequency of oscillation as the frequency of the 
forcing function. In fact, the frequency of the system response depends 
not only upon the frequency of forcing function but also upon the amplitude 
of the input. The principle of superposition no longer holds and the linear 
circuit characteristic is not applicable. In general, it can be said that 
the nonlinear system can have wide variety of almost periodic oscillations, 
the frequencies of which differ from the frequency of the applied force and 
vary with time as well as with different initial conditions. 

Sub-harmonic oscillation always existed with a frequency higher than 
the frequency of free oscillation with a particular amplitude of error signal. 
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The response curve which is a part of an ellipse or a hyperbola in the 

plane is always nearly parallel to the right side of the free oscil- 
lation response curve. A typical sketch of a soft saturation nonlinearity 
is shown in Fig. 4-1. 




Fig. 4-1 Fundamental and Sub-harmonic Response 

Curves for a Soft Saturation Nonlinearity 

Under some suitable conditions, there may exist any order of sub- 
harmonic oscillation, or only one type of order can exist depending on 
the characteristic of the nonlinearity and the input forcing function. 

Sub harmonic oscillation and response of nonlinear vibratory systems 
are often discussed in the field of mechanical system There is no 

body of knowledge applied to feedback control sexrvome chan ism. For the feed- 
back control system, there is a standard differential equation (3-6a), it 
is very similar to the equation for the mechanical system, the only dif- 
ference is that the forcing function of a feedback control system is a func- 
tion of input frequency. 

Although any order of sub-harmonic oscillation can exist in a non- 
linear feedback control system, in this paper only the sub-harmonic of order 
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2 and 3 are discussed* 



4-2 Sub-harmonic oscillation of Order 2 without Damping: 

Sub-harmonic oscillation of order 2 can exist in most nonlinear feed- 
back controls. When it occurs, the frequency of output (Co) or error 
(e ) is only 1/2 the frequency of input { 6Jr ) • A typical diagram from 
the computer experimental work shown on the "Input vs Output" plane and the 
waveforms of "Input and Output" are shown in Fig* 4-2 and 4-3, respectively* 
Recall the basic equation for forced oscillation in nonlinear system 
with un-damping Eq. (3-11). 

E + Kf(E) = -A Cos( tv. t + 9) (3-11) 

For the same type of restoring force function with chapter 3, leads: 

E + Ka^E + Ka^E^ = - A Cos{iv t + G) (3-12) 

Assume one solution of equation (3-12) is; 

E = E 2 C 0 S ^u-'t + E^ Cos t (^“1) 

Where amplitude of sub-harmonic of order 2, E^^ is the amplitude 

of fundamental frequency, since: 



" = 4 



' 1/2 



Cos yiOt - 



iV’ Cos cO t 



(4-2) 






t 4 . U=y^ r ^ ty, r ^ c T jf 



i ■ C< jr. r (4-3) 

We are only interested to the terms of second sub-harmonic and fund- 
amental frequency, inserting equations (4-1) to (4-3) into equation (3-12) 
the result is: 

(- r/y Ka^Ly L-y i- t (~ lO"B, i ho.h. 



•5 - . V I 

r r ^ Y , 1-y "o- 



. -.wr 



A’ >• . 



(4-4) 
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for a Sjsle:,!: f(2) --^2E - 0 .004of3 
K - .^ -0,25 

ei.(t):--r 2.50Cos( M. t - 0) 
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As discussed In Chapter 3, Section 3-3, if the frequency is larger 
then the frequency of free oscillation, the phase angle 9 is zero on the 
’’Error vs Frequency” plane, hence equation (4-4) leads: 

ilj- ^ (4-5) 



W' 



a - p 



(4-6) 



Equation (4-5) is the condition for existence of the sub-harmonic of order 



2. Solve 



j ^ _ _i iii ^ / yC 

~ ( j/(a) ^ ^ 1^1 J 

Since E^/2 real; which leads: 



(4-7) 



(X) < (4 h'd, t- V 

For the first iteration method, consider equal to zero: 

Equations (4-5) and (4-6) leads: 



(4-8) 



(4-9) 



t, ^ ft (4-10) 

Inserting equations (4-9) and (4-10) into equation (4-8); 

CD JZCCc'.-, t fi ic'd] < 0 (4-11) 

From Equation ('--5), it represents an ellipse or a phperbola in the E|^2 
vs ( 4 ) plane depending on the sign of a a^, also has a maximum when 
is less than zero, and a minimum when a^ is greater than zero. Actually, 
the suh-harmonic oscillation only existed within an interval of frequency. 
That means only one part of ellipse or hyperbola can exist, as shown in Fig. 



(4-4). 

4-3 Sub-harmonic Oscillation of Order 2 with Damping: 

Recall the basic equation (3-41) for a damped system: 

^ 3 

E + ^^1^ ^ ^^3^ “ ” ACos(o; t 4- 0) - BSin((t t + 9) (3-41) 
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Assume one solution of equatioii (3-41) is: 



t — £/ 



'■ ^ ic c- ' C/f.-; t 



(4-12) 



Since : 






^ ^ - - !' 3-f- ?. -4- "T t- *,.. , I* ^ 



(4-13) 

(4-14) 



^ 'A ' -^ 



t 4 - •'- 

f 'iS r ^ ^/A ^ ^ 



(4-15) 



Insertinc^ equations (4-12) to (4-15) Into equation ('^-41) and equal- 
ing the coefficients of Cos ^(A^t , Cos:t'^, and give a result of 

^ ^ (4-16) 

' iV r t >■ ^ >'-,t ^ 

A:ZAZ. Aj '\y^ uP (4-17) 

~t,jj(xr - v.n^yf lO t Ku> t,i) ^ (-^ £;'/, , 

- A a) ~ /3 c <>“■' (P (^" 1^) 



From equation (4-18), it is the coefficient of the term of Cos y c, 
where ^i/2 
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2 



(4-19) 



Solve for 

61/2 - ■ ' , _ ' (4-20) 

The value of ^^^{2 real therefore, the existence of 2nd order 

sub-hamnonic with damping: 

iij C J -f ^ L- k :* 'T . { \ ®3 * ^ (4-21 ) 

From equations (4-16) to (4-18) for the first Iteration method, for 
a^ equal to zero: hence: 

0^ s 4Ka^ (4-22) 

(ICa^i - SB - ^ FCos0 - FSinO (4-23) 

(Ka^ - oJ )E^^ ® FSinG - FCosG (4-24) 

Solve and E^^ from equations of (4-22) to (4-24): 



E 



lA 







-• V j 7 



'IB 



tL 



?f *■ 

/ '^>1 



t-f 



(4-25) 

(4-26) 



Inserting equations (4-25) and (4-26) into equation (4-21). A 
general equation for existence of 2nd order sub-harmonic oscillation 
with damping: 



0) :2. / 



( Z < 

1 lX> 7^ 't ^ K<-'( 



f Ka, h\ (' ^ 



/W ■* - -4 jL 



V , 
, ; : 



tbKcliY^l 



- ^ ^ ^ , c< . r f ^ 4 ^ J J ' , 






t 4 J • 



(4-27) 



For the case of G equal to 0 or n, i.e., SinG equal to zero equation 
(4-27) becomes: 




/V 






z 



'^4 ^ ■ 



'■ C k<^-\ 




(4-28) 
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From equations (4-11) and (4-28) they were shown that either the 
system with damping or without It the second order sub-harmonic oscilla- 
tion only existed with a frequency less than the frequency of two times 
the natural frequency of linear system for the case of less than zero, 
and greater than this frequency for the case of greater than zero, 

4-4 Sub-hannonlc Oscillations of Order 3 without Damping: 

Typical waveforms with jump phenomena of 3rd order sub-harmonic 
forced oscillation is shown In Fig, 4-5; (a) shows the transient from 
the 2nd order to the 3rd order sub-harmonic; (b) is shown reversing. 

For investigating the existence conditions of 3rd sub-harmonic forc- 
ed oscillation recall the basic equation of (3-11); and leads: 



E -f Ka^E + ^^^ 3 ^ * - A Cos(wt + 0) (3-12) 

Assume one solution of equation (3-12) is: 

E * Cos j iv t 4 E^ Cos^vt (4-29) 

in which is the amplitude of 3rd order sub-harmonic osclllationg 

hence: 

E » - ^ ^ 1/3 3 ^ ^ Cos 4v t (4-30) 

3 3 

' ' ^ * ®l/3^ " '^°*3 ■= 



13 2 

^^ 1/3 ^1 3 E j ) Coscvt 4*, ^,4,, 

From equation (4-29) to (4-31) and equation (3-12) leads: 



(4-31) 






'■ ^y. *■ '7-' 



(4-32) 



« A'i/ - *■ •• ^ 4 - ^ ' 

Equation (4-32) Is not equal to zero, hence; 



) ■ I f 



‘■-y ^ t,.. 

/ ■> 



/ 



('*-33) 



(4-34) 



Solve from equation (3-34). 
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i 



f 







(4*35) 



The value of be real, then: 

^ /v«. r c (4-36) 

Equation (4-36) is the condition for the existenceof sub-harroonic 
oscillation of order of q without damping system: 

From equations (4-32) and (4-33), by the first Iteration method, 



for a^ « 0; then: 



CC' 



9 Ka, 



(Ka^ - u )E^ * (o F 

Solve from equation (4-37) and (4-38): 

E s - F 

Inserting equation (4-39) in equation (4-36) the condition for a 3rd 
sub-hariTionir oscillation becomes: 



(4-37) 

(4-38) 

(4-39) 



0.} 



? •/ -V ~ * jjr tea if- ' 



Also the response equation of equation becomes: 









(4-40) 



(4-41) 



The equation of (4-41) represents an ellipse or a hyperbola in 
the vs a. plane depending on the sign of a^, also has a maximum 

when 0, and a minimum when A sketch of a 3rd sub- 

harmonic forced oscillation is shown in Fig, 4-4, 

4-5 Sub-harmonic Oscillation of Order 3 with Damping: 

Recall equation (3-41): 

E + E + Ka^E + Ka^E^ * - ACos(.ot + 0) - 3Sln(iv t f 0) (3-41) 

Assume one solution of equation (3-41) is: 



E * E* * ^Cos *Sj t 4* E-.CoS(X^t + E.^Sinw t 
1/3 3 lA IB 



(4-42) 
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Since : 



• 1 

E - i t - ^ Sin-t + £-_*- Cos t 

1/J ^ 



' ■ ■ ?i/3 ‘=<'' i '■ 



E,,.v Cos • t - - E,„ Sin- t 
lA IB 



(4-4-) 

(4-44) 



- I ^/3«l/3' * h//-lA * ^hA * 2E,3^Cos i..,t . 



4^1/3 ^ ^ ^^1/3 ^lA ^^lA ■*■ 

3'=1a''iB^> * 7 'iB <"lB^ * 2 E,/3' * Ej/) Sln» c * ..(4-45) 

4 

From equation (4-42) to (4-44) and (3-41), by equaling the coefficients 
of Cos ■“ % t. Sin t. Cost t, leads: 






' t ■ •- > - -.u > 



r^i> ■ ^: ) 

(A -2/ 

Rearrange equation (4-46) and solve it for terms of 



.fr.' V (t:, 



/ « ) r^/i 

~ // --4 -r ;. 

- A ; t h 



I ♦ • 



(4-46) 

(4-47) 

(4-48) 



Eia, and Ejg: 






^ . h :> 









here should be renl^ hence: 






or: 



i ^ ^1/ >7 



u 



y r .• (c" 

( . ^ r ..-r/ 

Equation (4-51) is the condition of 3rd sub-harmonic oscillation 
existing in saturation nonlinearity with damping system. 

By the iteration method, equation (4-46) to (4-48) leads! 

.1 - 9Ka, 



* V ; 



r 



(4-49) 



(4-50) 



(4-5i) 



(4-53) 



(Ka^ 


- cv )Ej^^ + UX Ejg = - 


1 FCos 9 - ^ 


'FSine 


(4-54) 


(Ka^ 


-A )E^^ - y- = , 


FSin 0 - V / 


FCosG 


(4-55) 
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Solve equation (4-53) to (4-55), 



E ~ 
lA J 



(4-50) 



‘IB 



- . • ■ (4-57) 

I •* f 

Substituting equations (4-56) (4-57) Into equation (4-52) and be- 
comes in a general form: 






, L i / 






f‘ t 

i ^ r 



l- J 7 f <• i ' r y. 



-r -V r 



(4-58) 



In the case of 6 equal to 0 or rt, l.e, SinG equal to zero equation (4-58) 
becomes : 



A ^ 3 




IJ-fXhU 




A 






/ 



/ ^ 






(4-59) 



In conclusion, the condition for existence of sub-harmonic oscilla- 
tion in a nonlinear feedback control system depends on the characterlcs 
of the nonlinear element and the amplitude of input. The response curve 
of any order sub-harmonic oscillation with a saturation nonlinearity is 
either an ellipse or a hyperbola which depends on the characteristic of 
nonlinearity. Also the frequency of a sub-harmonic oscillation existed 
may be either a maximum when a soft saturation nonlinearity is presented 
or a minimum when a hard saturation nonlinearity is presented. The limit 
value of minimum or maximum is almost equal to the number of sub-harmonic 
order times the natural frequency of linear system. 

4-6 Analog Computer Analysis for a 2nd Order Sub-harmonic Oscillation 
System: 

Consider a system with a block diagram shown in Fig. 4-6 t 
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Fig. 4-6 Block Diagram of a 2nd Sub-harmonic 
Oscillation System. 

Assume the restoring force function of nonlinearity is: 

f(E) = 2 e - 0.0048E^ (4-59) 

The characteristic curve for equation (4-59) is shoxvn in Fig. 4-7, in 
which: 

E * 1.00; 2.50; 5.00; 7.50; 10.00; 12.50; 15,00. 

f(E) = 1.995; 4.295; 9.400; 12.97; 15.20; 15.40; 14.80. 

For the simulation of restoring force function of nonlinearity from 
computer, it is approximated from two straight lines, one of which with 
a slope of 2; and the other with a slope of zero. The circuit setup in the 
computer is identical with the circuit shown in Fig. 3-15. The only dif- 
ference is the value of bias voltage to be used. The characteristic curve 
of restoring force function from the computer is shown in Fig. 4-8. 

An operational diagram setup in the analog computer is almost the 
same with the diagram of Fig. 3-20, only one difference is that the value 
of constant to be used. 

The results from the computer are condensed in Table 4-1 and Fig. 4-9. 
Only the response curve is of interest, hence, only one response curve in 
the error vs frequency plane is presented. It is shown that there are two 
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jumps, one is the fundamei»tal frequency Jump, and the other is the second 
order sub-harmonic frequency jump. These jump phenomenas are shown sche- 
matically in the diagram of Fig. 4-9, 



Table 4-1: 



Data for a 2nd Order Sub-harmonic Oscillation Response 
Curve, in which; f(E) * 2E - O.OOASE^; K » 4^</« 0.25; 
and 9^(t) « 2,5 Cos( - t *f 9) 



(a) Frequency Increasing: 



quency of 


Amplitude of 


Amplitude o 


Input 


1st Harmonic 


2tid Sub-har 


( ) 


( e,) 


( ^ ) 


0.628 


0.22 


/ 


0.754 


0.33 


/ 


0.880 


0.45 


/ 


1.000 


0.50 


/ 


1.130 


0.63 


/ 


1.260 


1.00 


/ 


1.380 


1.15 


/ 


1.510 


1.50 


/ 


1.640 


2.00 


/ 


1.760 


2.50 


/ 


1.890 


3.00 


/ 


1.950 


3,25 


/ 


1.980* 


10.50 


/ 


2.010 


10.00 


/ 


2.200 


9.00 


/ 


2.510 


7.30 


/ 


2.810 


6.00 


/ 


2.950 


/ 


/ 


3.140 


5,00 


/ 


3.460 


4.50 


/ 


3.770 


4,00 


/ 


4.090 


3.75 


/ 


4.400 


3.50 


/ 


4.710 


3.25 


/ 


5.020 


/ 


4.75* 


5.150 


/ 


4.50 


5.270 


/ 


3.50 


5.340 


2.50 


/ 



Rema rk s 



^*lst harmonic 
upper jump 



''2nd sub-harmonic 
upper jump 
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(b) Frequency Decreasing: 



Frequency of 
Input 



Amplitude of 
1st Harmonic 



( - ) 


( - ) 


( X- 


5.340 


2.50 


/ 


5.270 


/ 


3.50 


3.150 


/ 


4.50 


5.020 


/ 


4.75 


4.710 


/ 


5.50 


4.400 


/ 


6.50 


4.000 


/ 


7.80 


3.770 


/ 


9.50 


3.460 


/ 


11.00 


3.140 


'/ 


13.50 


2.950 


/ 


14.50* 


2. 810 


6.00 


/ 


2.510 


7.25 


/ 


2.200 


8.75 


/ 


2.010 


10.00 


/ 


1.890 


11.00 


/ 


1.760 


12.25 


/ 


1.640 


13.00 


/ 


1.510 


15.00 


/ 


1.380 


17.00 


/ 


1.260 


21.20 


! 


1.130 


24.00 


/ 


1.000 


25.50 


/ 


0.880 


0.45* 


/ 


0.754 


0.33 


/ 


0.628 


0,22 


/ 



Amplitude of 

2nd Sub-hnniionlc 
) 



Remarks 



^2nd s ub - he mion i ' 
lower jump 



-J^lst Harrnonic 
lower jump 



79 



I 



I 



I 

I 

( 

I 

I 

I 

I 

I 

I 

I 



i 



I 

i 




80 



I 

I 

I 



I 

I 



I 

I 



I 

I 




4*7 Analog Computer Analysis for a Ird Order Sub-harrnonlc Oscillation 
System: 

Consider a systen. with a block diagram shown in Fig* 4-10: 



in which assume the nonlinearity of system is an "ON-OFF” element, a 
sketch of the characteristic for the computer simulation is shown in 
Fig, 4-11, it is appro/cimated from three straight lines; one of which 
with a slope of 15,0; one of which with a slope of 0,1; and the other 
is zero: 

The circuit setting up in the computer is shown in Fig, 4-12, all 
components to be used are schematically indicated with the diagram. The 
characteristic curve from the computer is shov;n in Fig, 4-13, 




j / ^ 



r 



Fig. 4-10 Block Diagram of a 3rd Sub-harmonic 
Oscillation System 



H- 







r 



Fig, 4-12 Function Generator Diagram for an 
"ON-OFF" Nonlinearity 
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An operational diagram tor setting up the analog computer Is shown 
in Fig* 4-14, in which the value of all components to be used are 
schematically Indicated in the diagram. The scaling factor for the 
computer setup is the same as diagram of Fig. 3-20 to be used. 

The results from the computer are condensed in Table 4-2, and with 
a response curve shown in Fig. 4-15. It is shown that, there is a pair 
of jumps for each harmonic oscillation response curve, these jump pheno- 
menas within frequency range is accompanied shown with the diagram of 
Fig. 4-15. 
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Treble 4-2: Data for a 3rd Order Sub-harmonic Oscillation Response 

Curve from an *’ON-OFF** Nonlinearity System and; ^ » 0.4; 
K * 10.0; 0v^(t) XT 8.0 Cos( ot 4 0) 

(a) Frequency Increasing: 



Frequency of 


Amplitude of 


Amplltdue of 


Amplitude of 


Bema rks 


Input 


1st Harmonic 


2nd Sub-har. 


3rd Sub-har, 




(CO ) 


( fc/ ) 


( 


( 




0.628 


0.60 


/ 


/ 




0.880 


1.25 


/ 


/ 




1.130 


2.25 


/ 


/ 




1.380 


3.50 


/ 


/ 




1.640 


4.75 


/ 


/ 




1.760 


5.75 


/ 


/ 




1.890 


7.50 


/ 


/ 




2.010* 


50.00 


/ 


/ 


*lst Harmonic 


2.140 


46.00 


/ 


/ 


Upper Jump 


'>.270 


42.50 


/ 


/ 




2.390 


39.00 


/ 


/ 




2.510 


35.00 


/ 


/ 




2.630 


31.50 


/ 


/ 




2.760 


28.00 


/ 


/ 




2.890 


26.00 


/ 


/ 




3.010 


22.50 


/ 


/ 




3.140 


20.00 


/ 


/ 




3.260 


/* 


/ 


/ 


^Transient and 


4.270 


/* 


/ 


/ 


unstable 


4.390“ 


/ 


/ 


/ 


^Transient trom 


4.520 


/ 


/ 


/ 


1st Harmonic 


4.650 


/ 


/ 


/ 


to 2nd Sub-har. 


4.770 


/ 


/ 


/ 




4.990“ 


/ 


/ 


/ 




5.020* 


/ 


/ 


/ 




5.140 


/ 


20.00 


/ 


^2nd Sub-ha i. 


5.270 


/ 


18.00 


/ 


Upper Jump 


5.400 


/ 


16.70 


/ 




5.520 


/ 


16.00 


/ 




5.650 


/ 


15.00 


/ 




5.770 


/ 


13.80 


/ 




5.900 


/ 


13.50 


/ 




6.030 


/ 


12.70 


/ 




6.150 


/ 


12.30 


/ 




6.280 


/ 


11.90 


/ 




6.400 


/ 


11.30 


/ 




6.530 


/ 


11.00 


/ 




6.650 


/ 


10.60 


/ 




6.780 


/ 


10.20 


/ 




6.910* 


/ 


10.00 


/ 


*3rd Sub-h r* 


7.030 


/ 


/ 


35.00 


Upper Jump 


7.160 


/ 


/ 


32.50 




7.280 


/ 


/ 


28.50 




7.410 


/ 


/ 


27.20 




7.530 


/ 


/ 


26.00 
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(b) Frequency Decreasing: 



Frequency of 
Input 
) 

7.530 
7.410 

7.280 
7.160 

7.030 
6.910 
6.780 
6.650 

6.530 

6.400 

6.280 
6 . i 50 

6.030 

5.900 

5.770 
5 . 650 * 

5.520 

5.400 

5.270 

5.140 
5.020 

4.900 

4.770 
. 4.650 

4.520 

4.390 
4 . 27 C 

4.140 
3 . 260 ' 

3.140 

3.010 

2.890 

2.760 
2.630 

2.510 

2.390 

2.270 

2.140 

2.010 

1.890 

1.760 
1 . 640 * 

1.510 
1.380 
1.130 
0.880 
0.628 



Amplitude of 
1st Harmonic 
(4 ) 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

20.00 

22.50 

26.00 

27.50 

31.00 

35.00 

39.00 

42.00 

45.50 

51.00 

55.00 

52.50 
5.00 

4.25 
3.50 
2.30 

1.25 
0.60 



Amplitude of 
2nd Sub-har. 



/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

13.75 

15.00 

16.00 
16.70 

17.50 

19.50 

21.50 

23.00 

25.00 

32.00 

36.00 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 



Amplitude of 
3rd Sub-har. 

26.00 

27.20 

28.50 

30.50 

32.50 

35.00 

36.50 

37.50 

40.00 

42.50 

43.75 

47.00 

50.00 

52.50 

56.00 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 



Pema rks 



*3rd Sub-har. 
Lower Jump 



®*2nd Sub-har. 
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CHAITEK V 



FORCED OSCILI/.TION FOR THE HIGHER ORDER 
FEEDBACK CONTROL SYSTEM 
5-1 General Description ; 

Up to the Chapter 4, only forced oscillations in the case of the 2nd 
order system are discussed. In the extension to higher order systems, 
consider a block diagram shown in Fig. 5-1; 



o I 






C/ A/ 






) 






in which G(s) 
of the system. 



Fig. 5*^1 Block Diagram of a Feedback 
Control System 

is an equivalent transfer function of the linear element 
it can be written in tVie form: 



G(s) 



K 



F(s) 



( 5 * 1 ) 



and with a result for the system; 

P(s) E + Q(s) Kf(£) = P(s) (5-2) 

Where I(s) may be any order of s; as: 

P(s) ss s” + h.s” ^ +■ + b (5*3) 

12 n 

Slmllarlly Q(s) may be 

_ . , m m-1 m-2 , . 

Q(s) s s + Cj^s + C 2 S + -f (5*4) 

In general the frequency responsive element 3(s) following the non- 
linearity is of such a nature as to attenuate the high frequency^, that 
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me^rns the order of P(s) should bu higher than the order of Q(s), The 
order of system is n. Inserting equations (5*3) and (5-4) into (5-2) and: 



■dl 

ic. ' ‘ ./r 




(!>■’ ' 

n't ' ’■ 






,n 

('u ^ ^ ^c>- I '■ ’ ) (5-5) 

Equation (5-5) is a general differential equation form of the 

forced oscillation of any order feedback control system. 

5-2 Harmonic Linearization and Iteration Method for Investigating the 
Equivalent Gain of a Nonlinearity: 

After a restoring force function of the nonlinear element has been 

defined, the equivalent complex gain, N(Ej) of the nonlinear element is 

determined by the ratio of the complex amplitude of the first harmonic 

of the output to the amplitude of the harmonically varying input quantity. 

Recall the block diagram of Fig. 5-1, if the system is a 2nd order 

3 

system and with a restoring forcing function; f(E) s a^E + a^E , and: 



K N(E) E s= 0 
s (s ^ 

or: 

E + E + K N(E) E = eV + 
Assume the input function is 

li’, ss F Cos (a t + 0) 
and one solution of E: 



E = Ej^ Cos (V t 



and with a result: 

( >:" E^ + K N(Ej)Ej)Cosco t E^Slnu t 



-A Cos(i;, t + 9) 



(5-6) 



(5-7) 



(5-8) 



(5-9) 



-BSin(£^t + 0) (5-10) 

Comparing the result with the equation of (3-42) the equivalent gain: 



N(Ej) 




(5-11) 
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For a general lorm expression; the equivalent gain of a non- 
linearity In terms of the airplitude of Input Is devtl« ped by Ya, Z, 

Tsypkin, The result for a syrnmetrlcal nonlinearity with input of equation 
(5-9) is: 

N(Ej) = [ f(Ep + f(Ei/2)/ (5-12) 

If more accuracy is required, the equivalent gain N(E^) may be; 

N(Ei> - [ f(Ep 4 f(E^/2) 4- ,3 f( v3 £^/2), (5-13) 

For the case of a nonlinearity with a restoring force function^ 

3 

f(E) « a^E + a^E , the results from equations (5-12) and (5-13) are ex- 
actly the same as equation (5-11) • Hence, for any nonlinearity the case 
of first harmonic input the equivalent gain can be expressed either in the 
form of equation (5-12) or (5-13). 

Equation (5-12) or (5-13) not only substantially simplifies the 
calculation for a second order forced oscillation, it is also probably 
a most important for investigating self oscillation or forced oscillations 
in a higher order system with single or multiple nonlinearities in a 
feedback control system. 

5-3 Forced Oscillations for a Higher Order Feedback Control Systems 

Consider a system with a linear transfer functions 

G(s) = (5-14) 

s(s + Pj)(s + P 2 )(s 4- p^) 



or in the form: 



G(s) 



K(s + c) 



4 3 2 

s 4- bjS 4- b 2 S 4- b^s 



(5-15) 



From the equation (5»5), the differential equation for the system 



of equation (5-15): 
.4 



4 b -i- “77 ♦ Of(E) 

( ~ ^ 1 ~ •’2 ^2 * ^3 dr 

dt 



dt 



df 
= ( 



4 3 2 

±- ^ K il- ^ h £-9 ^ K -±~> 

.4 ^1 . 3 ^2 dt ^ ^3 dt 

dt dt 



(5-16) 



oi 



in which, 



E Ej^Cos tw L 

: = FCos(^t + 0) 



and ; 



1 3 

f(E) « E^Cos ^ ^3^1 (^Cosvot + Cos3 v t) 

Inserting these equations into equation (5-16) and with a result! 



-t: fr,ay t r ^ r<z,c c/ j ^ 

T ( h, L. CO — ^ /i .0/. 

r Y K a J ~ ^/iVj >' 

~ i-W I Ci/ { oi '' ~h^)cvii^ r [h,l>y'~ hj, ’i. Jc 



- /'H. / St) ( i,\) ' 0 ■■ [b/i^ ~ '' , ) Cc: Sy j f' v’ 

(5-17) 

Neglect the higher order harmonic terms and equate the coefficients 
of Cosset and Sina t: 

ififeo '' - h .Ci.' ' r A it, c ' i c ^ 

“ / : -t ( b.i‘ I (5-18) 

^>Cb. “(A-. .. r- i 

- t ->r <)c\) ~ i hii')'' (5-19) 



Squaring equations (5-18) and (5-19) and adding; 



Er<- 



it ' 



Ki'^ 



i •- 



t<ai L :r, b. /i.' b^) 7,f2, r ' 



— f tv / iv '■ ( I ' ' r T rbjti/ ii) J 

(5-20) 

Equation (5-20) is the equation of response curve of the system of 
equation (5-15). If the constants of system are know, a response curve 
in the s?, vi^ plane can be plotted. 

For the investigation of conditions for existence of forced oscil- 
lation rearrange the equation (5-20) and put : 
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I 

I 

I 

I 
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I 
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X 



= a. 




Hence : 



p ' . . / = A 






(5-21) 



or; 



V'* T 









- 




h 

t [ ^ - r, r ^ ‘ ^ - (5-22) 

The condition for existence of forced oscillation Is found from 
equation (5-22) by means that the solution of should be real. It Is 
to be noted that the value of X Is the equivalent gain of nonlinearity, 
its value should be larger than zero,, as It was proved in Chapter 3. 

For the phase Investigating recall the equation (5-19) and re- 
write : 

V * - h - ) u ^ • > 

hi* 

I. r. U-i ' - I !-■ (5-?3) 

or: 



C O -./ / 



- ' ' 'j. i : J -■ 

I [U I .. I- " 






hence : 



(9 






/ ' c t ^ ) r- I i , 






(5-25.) 



in which the value of; 

^ =. -an. - h , ) l< ^ )j (5-/.v) 

Equation (5-25) is shown that the value of 6 Is a function of .. 
and the amplitude of first harmonic of error signal, if the system 
constants is fixed. 
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If we use the equivalent gain for the system investigation, and 
with a result: 

4 — ^7 ( 5 . 27 ) 

/ t /V. t. ; H ; 

or : 

^ (5-28) 

Where is the amplitude of first harmonic of error signal of forced 
oscillation. 

Equation (5-28) is similar to a linear system, the techniques for 
linear system are applicable. 

5-4. Analog Computer Analysis for a Higher Order System: 

Consider a feedback control system with a transfer function and a 
nonlinearity are shown in Fig. 5-2: 





[ ^ 17 ^ 






> 0 


?) '-• >j -s...-. j/:--.'! 




rc \;6 ^ ' 








4 



Fig, 5-2 Block Diagram for a Typical 
Feedback Control System 

in which the restoring force function is: 

f(E) « 2E + 0.020E^ (5-29) 

The characteristic curve of restoring force function of nonlinear- 
ity from the calculation is shown in Fig. 5-3. Where: 

E = 1.00; 2.00; 3.00; 4.00; 5.00; 6.00; 7.00; 7.50; 

f(E) = 1.98; 3.84; 5.46; 6.72; 7.50; 7.68; 7.15; 6.60: 
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C'5-3 Character! stiq Cui’vc of 
f(E)=*2B ~ 0.0206-^ Fnst?! f 
Calovilation 
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For the sim»ilation oi the restoring force £»jnctlon of nonlinearity 



from the Computer, it is approximated by three straight lines, one of 
v^hich with a slope of 2; one of which with a slope of 1; and the other 
with a slope of zero. The circuit setup in the computer is shown in 
Fig. 5-4, all components used are schematically in the diagram. 

The characteristic curve of restoring force function from the com- 
puter is shown in Fig, 5-5; (a) shows the waveforms of the input and the 
output; (b) shows the characteristic curve of the nonlinearity. 



The operational block diagram and the operational circuit diagram 
are shown in Fig. 5-6 and Fig. 5-7. The all components used in the 
operational circuit diagram are sechemat Ically In the diagram. 

The results from the computer are condensed in Table 5»1, and with 
response curves shown in the plane and plane are sho^vn In Fig. 5»8 
and Fig. 5-9 respectively. 





t 0 




L 



Fig. 5-*4 Analog Computer Setup for 

f(E) s 2E - 0.020 e 3 Simulation 
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Fig. 5-6 Operational Block Diagram setup 
in the Analog Computer 
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Table 5-1* Data for Re^^ponse Cuv\e from Computer for a Hlghet Order 

System; in which f(E) - 2E - 0.020 e 3, 6/ (t) - 3.0 Cos(^.’t +0) 

(a) Frequency Increasing; 



Frequency of 


Amplitude of 


Phase of 


Amplitude 


of 


Phase of 


Input 


Error Signal 


Error 


Output 




Output 


(iO ) 


( ^/) 


( 4 ) 


( ) 




( 4 ) 


0.314 


0.40 


270.00 


3.10 




0.00 


0.377 


0.55 


270.00 


3.30 




0.00 


0.440 


0.75 


262.00 


3.40 




0.00 


0.503 


0.90 


255.00 


3.50 




5.80 


0.565 


1.10 


246.00 


3.60 




9.00 


0.628 


1.25 


241.00 


3.80 




12,50 


0.755 


1.70 


235.00 


4.30 




20.00 


0.880 


2.50 


220.00 


5.00 




24.50 


1.000 


7.40 


105.00 


6.30 




107.50 


1.130 


6.90 


64.00 


5.10 




120.00 ' 


1.260 


6.50 


37,50 


4.25 




137.50 


1.380 


6.10 


25.00 


3.50 




155.00 


1.510 


5.60 


15.00 


3.00 




165.00 


1.640 


5.25 


10.00 


2.50 




172.00 


1.760 


4.90 


7.50 


2.00 




180.00 


1.890 


4.50 


5.00 


1.70 




186.00 


2.200 


4.00 


0,00 


1.20 




195.00 


2.500 


3.50 


0.00 


0.85 




210.00 


2.810 


3.30 


0.00 


0.60 




220.00 


3.140 


3.20 


0,00 


0.40 




232.00 


3.460 


3. 10 


0.00 


0.30 




245.00 


3.770 


3.00 


0.00 


0.25 




260.00 


4.090 


3.00 


0.00 


0.20 




270.00 


(b) Frequency 


Decreasing 










Frequency jf 


Amplitude of 


phase of 


Amplitude 


of 


Phase of 


Input 


Error Signal 


Error 


Output 




Outr-ut 


( a-' ) 


( <-)) 


( 4 ) 


( ) 




( S ) 


2.510 


3.50 


0.00 


0.78 




210.00 


2.200 


4.00 


0.00 


1.20 




195.00 


1.890 


4.50 


5.00 


1.65 




185.00 


1.760 


4.90 


7.50 


2.00 




180.00 


1.640 


5.25 


10.00 


2.50 




171.00 


1.510 


5.70 


15.00 


3.00 




162.00 


1.380 


6.10 


25.00 


3,50 




155.00 


1.260 


6.60 


37.50 


4.25 




137.00 


1.130 


6.60 


64.00 


5.00 




121.00 


1.000 


7.40 


105.00 


6.25 




106,00 


0.943 


/ 


/ 


6.80 




78.00 


0.880 


6.70 


150,00 


6.00 




25.00 


0.754 


1.80 


235.00 


4.30 




20,00 


0.628 


1.25 


248.00 


3.80 




16.00 
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CHAPTEF VI 



CONCLUSIONS 

The aim of this dissertation was to present a method for investigat- 
ing the conditions for existence of the forced oscillations of a nonlinear 
feedback control system. This method is called "harmonic linearization 
method", it is based on the harmonic balance and Iteration method. From 
this method, the conditions for existence of the fundamental frequency 
and the sub-harmonic of order 2 and 3 are investigated, and also a equi- 
valent gain of the nonlinearity in terms of input amplitude and frequency 
has been developed. 

The term of the "response curve" as used here is defined as the combina- 
tion of two components, one is the attenuation frequency; It is either 
shown in the "Error vs Frequency" plane or "Output vs Frequency" plane, 
and the other is the phase frequency response, it is with a difference 
value from the different representation of attenuation frequency response. 

The results of this analysis show that a nonlinear feedback control 
system may possess more than one type of forced oscillations, which type 
will exist, it is depending on the characteristic of the system and the 
characteristic of input forcing function. 

For any nonlinear feedback control system, the response curve of 
either fundamental or sub-harmonic forced oscillations, it is usually a 
sudden jump in amplitude and phase when a jump transient take place. The 
range of jump frequency depends on the value of damping of the system and 
the amplitude of the input; usually it is increased as the damping decreas- 
ing, and decreased as the input amplitude decreasing, A vertical tangents 
locus investigating for calculating the jump frequencies in terms of system 
constants and the amplitude of input was investigated. 
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Due to the resonance phenomena and the linear approximations of re- 
storing force functloti, there is some difference between the theoretical 
and the experimental response curves. It ir only exLstant in the jump 
frequency range, other than this range of frequencies, it is very close 
to each other. 

Sub-harmonic oscillations are undesired for a nonlinear feedback 
control system, it is actually an unstable phenomenon. Under suitable 
conditions, a multiple order sub-harmonic forced oscillation may exist 
in a system with a constant input forcing function. There is also a 
jump phenomena for each order of sub-harmonic oscillations. The condi- 
tions for the existence of the 2nd and 3rd order sub-harmonic forced 
oscillations are investigated, and also theoretical response curve equa- 
tions for each case are developed. 

It is possible an important advantage that the design and analysis 
of a nonlinear feedback control system by using the harmonic lineariza- 
tion method, by means for Investigating an equivalent gain of the non- 
linearity. In this case, all principles and techniques applied for linear 
system are applicable. 
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